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Abstract 

In this paper, we calculate the explicit expression for the fugacity for three-dimensional 
hard-sphere Bose and Fermi gases from their equations of state in isochoric and isobaric 
processes, respectively, based on the mathematical result of the boundary problem of an 
analytic function — the homogeneous Riemann-Hilbert problem. 

1 Introduction 

The physical problem. The equation of state for a quantum gas can be obtained by eliminating 
the fugacity z between the equations 

n=-^^^lnH(r,z,y), (2) 

where P is the pressure, InH the grand potential, n the particle number density, and V the 
volume. Concretely, for example, in isochoric processes, solving z = z{n,T) from eq. ^ and 
substituting it into eq. ([U give the equation of state 



(3) 

z=z{n,T) 



and in isobaric processes, solving z = z (P, T) from eq. ([T]) and substituting it into eq. ([2j) give 
the equation of state 



(4) 

z=z{P,T) 



However, only in some simple cases, such as classical ideal gases [T], two-dimensional ideal 
quantum gases, etc., the explicit expression of the fugacity z can be obtained by solving eq. 
dll) or ([U). In most cases, the fugacity can only be obtained approximately, e.g., at high- 
temperature and low-density or low-temperature and high-density limit. In this paper, we will 
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calculate the explicit expression for the fugacity for three-dimensional hard-sphere Bose and 
Fermi gases by directly solving eqs. ([2]) and ([I]) for isochoric and isobaric processes, respectively. 

The hard-sphere g simplified model, is of great value for investigating the more 

general theory of interacting gases and can be extended to some more general cases. This is 
because a particle that is spread out in space sees only an averaged effect of the potential and, 
thus, often a complete knowledge of the detailed interaction potential is not necessary for a 
satisfactory description [2]. Or, from the viewpoint of quantum mechanics, due to the low col- 
lision energy of collisions among the gas molecules, the shape-independent s-wave contribution 
dominates. The equations of state for three-dimensional hard-sphere Bose and Fermi gases are 
presented in Ref. [3] by the binary collision expansion method. 

The mathematical method. The key idea of the method for solving the explicit expression 
for the fugacity from eqs. ([2]) and ([1]) is based on the boundary problem of an analytic function 
— the homogeneous Riemann-Hilbert problem [3]. In a Riemann-Hilbert problem, one seeks 
to find a sectionally analytic function ip under the boundary condition 

^^+(0 = ^(0^- (0+5(0, Cei, (5) 

where L is a union of a finite number of smooth simple arcs, (C) and (C) are boundary val- 
ues of (0 on the left and right of L, G {() and g (C) are functions satisfying \f (C2) — / (Ci)l ^ 
A 1^2 — CiT with A, fi > 0, the Holder condition, and G (0 7^ everywhere on L. The func- 
tion ip (0 being a sectionally analytic function means that ip (0 is analytic in each region not 
containing points of the boundary L and is continuous on L from the left and from the right, 
excepting possibly some ends of L and near such ends the function (0 should satisfy 

constant 

where is a constant corresponding to the A;-th end Ck- The homogeneous Riemann-Hilbert 
problem is a Riemann-Hilbert problem with g (0 =0. In a homogeneous Riemann-Hilbert 
problem, the problem is converted into finding a function -0 (0 from the jump on the two sides 
of the boundary L,G{C)=i'^ (0 /V'" (0 • 

Take an isochoric process as an example. In this case, the fugacity is determined by eq. ([2]). 
The first step is to introduce a real function whose zero gives eq. ([2]). Analytically continuing 
this real function to the complex plane gives a complex function, and, clearly, the zero of the 
real function becomes some of zeros of this complex function on the real axis. Therefore, the 
fugacity z becomes a real zero of this complex function, and the problem of solving the fugacity 
is converted into a problem of seeking the real zero of such a complex function. The key step 
toward solving the zero relies on the homogeneous Riemann-Hilbert problem. Based on the 
Riemann-Hilbert problem, Leonard calculates the explicit expression for the fugacity for ideal 
Bose and Fermi gases [5|. More details about this method are described in Sec. [2j 

A description of the method is given in Sec. [2l In Sees. OandSl we calculate the explicit 
expression for the fugacity for hard-sphere Bose and Fermi gases. Sec. [3] for isochoric processes 
and Sec. S]for isobaric processes. The conclusions are summarized in Sec. [5l 

2 The method 

In this section, we give a description of the method. 



IV' (01 < 1^ afc < 1, (6) 
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2.1 The formal explicit expression of the fugacity 

Take an isochoric process as an example for describing the method. 

In an isochoric case, what we want to do is to solve an explicit expression for the fugacity 
z from eq. ([2]). Introduce a real function 

* (^) = -^^Yz ^iT,z,V)-l^f {z) - 1, (7) 

where f {z) = [l/nV) z{d/dz)\n'E{T,z,V). Clearly, eq. ([2]) corresponds to ^ {z) = 0. In 
other words, the fugacity z is a zero of the function ^'(-z), and the problem of solving the 
fugacity z is converted into the problem of seeking the zero of ^ {z). 

For seeking the zero of ^ (z), we first analytically continue the function ^ [z) to the whole 
complex ^-plane. The analytically continued function is 

^(C) = /(C)-i, CgC. (8) 

The fugacity z that we want to solve is a zero of the real function ^ [z) and is, of course, 
a zero of the complex function ^ [Q) on the real axis. Usually, the function ^ {(,) has not 
only one zero. Besides the zero corresponding to the fugacity z, there are still other zeros in 
the complex ("-plane, denoted as Wj, (i = 1, 2, • • • n^; — 1), where is the total number of the 
zeros. Moreover, ^ {Q) may also has singularities in the complex (^-plane. In the following, 
we will show that in our case, ^ {Q has no isolated singularities, and all its singularities are 
non-isolated singularities, forming some arcs in the complex (^-plane. Such arcs, which form a 
boundary of the analytic region of ^' (C); will be denoted hy Li, {i = 1,2, - ■ ■ ,mL), where rriL 
is the total number of the arcs. That is to say, ^ (C) is analytic in each region not containing 
points of the arcs Lj and has n^^ zeros; at least one of the zeros of ^ (() is on the real axis. 
Therefore, we can express ^ (C) in the following form: 

Uui — l 

*(C)=r?(C-^) n (C-^.)^(C), (9) 

i=l 

where is a constant and $ (C) is a function vanishing nowhere in the complex (-plane and 
analytic in each region not containing points of the arcs Lj. In principle, except the fugacity 
z, if we know all other — 1 zeros of $ ((), uji, the explicit expression of z can be formally 
expressed from eq. ([9]), 



c 



1 



(10) 



Then, the problem of solving z is converted into the problems of solving the function $ (() and 
finding the zeros Wj. 

2.2 $ {() and the fundamental solution of the homogeneous Riemann-Hilbert 
problem 

The function ^ (() can be determined with the help of the homogeneous Riemann-Hilbert 
problem. 



3 



In a homogeneous Riemann-Hilbert problem, if the jump of a function on the boundary is 
known, we can determine the function up to an arbitrary polynomial. For seeking the function 
we first need to solve the fundamental solution of the homogeneous Riemann-Hilbert 
problem, denoted as v?(C)- The fundamental solution ip{C) is such a solution that (C) and 
its reciprocal l/<p{Q both are sectionally analytic functions. More concretely, a fundamental 
solution has no zeros and isolated singularities; all its singularities lie on some arcs, forming a 
boundary for the analytic region; at each end of the arcs, c^, the degree of divergence of the 
solution (p{C) and its reciprocal is less than 1, i.e., |v'(C)| < constant / {l!^ — Ck\'^'' and 

\l/(p (C)| < constant/ \C, — Cfc|"*, where ak and are constants less than 1. The fundamental 
solution, generally speaking, is not completely determined by the above conditions; they are 
divided into some classes according to their behaviors near the ends, and any class of which 
can be chosen as the fundamental solution of the homogeneous Riemann-Hilbert problem. In 
this paper, we choose the fundamental solution satisfying 

, constant „ ^ , 

|y(C)l < |._ la, . 0<afc <1. (11) 

Generally speaking, the function <I> ((^) that we want to find is not a fundamental solution 
though it has no zeros and isolated singularities and all its singularities lie on the arcs, Li 
(i = 1, • • • , m/,), since the degree of divergence of $ (C) or 1/$ (C) at the ends of such arcs may 
not be less than 1. However, in terms of the fundamental solution ip the function <I> (C) can 
always be expressed as 

^{C) = ^r^ > (12) 

n - ^i^f' 

k=l 

where n is the number of the ends that are different from infinity and (3^ is a constant determined 
by both the degree of divergence of $ ((^) and the degree of divergence of the fundamental 
solution {() at the A:-th end c^. For a chosen fundamental solution if {() (in the present case, 
the fundamental solution is chosen to satisfy eq. (|lip ). f3k should ensure that the behaviors of 
the two sides of eq. (jl2p in the neighborhood of the k-th end of the arcs Li are the same. For the 
end at infinity, we need not pay special attention since at such an end the degree of divergence 
of the function encountered in the present case is less than 1. Concretely, near infinity, the 
behavior of the function ^ (C), defined by eq. lies on the asymptotic behavior of the Bose- 
Einstein integral or the Fermi-Dirac integral near infinity. Near the point of oo, the Cauchy 
principal value of the analytically continued Bose-Einstein integral go- (C) ~ — (InC)*^ / (o")] 
[8], and then ^(C)|(^_+oo ~ (l^C)'^) where A; is a constant; near the point of — oo, the Cauchy 
principal value of the analytically continued Fermi-Dirac integral (C) ~ [In i~C)V I V'^ ('^)]! 
and then ^ ^ \Lu{-Qf . 

The fundamental solution ((^) can be obtained based on the result of the homogeneous 
Riemann-Hilbert problem. From eqs. ([9|) and (jl2p . we can see that the jumps of the functions 
^ (C), ^ (C)) and Lp (C) at the two sides of the boundary, which consists of arcs Lj, are the same, 

<^+(C) $+(C) ^+(0 



^- (C) (C) ^- (C) 



G(C), CG^*,(^ = l,2,---mi). (13) 
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According to the homogeneous Riemann-Hilbert problem, from the jump G {Q), the fundamen- 
tal solution if (C) can be determined: 

n 

<^(C) = e^(^)n(C-Cfc)"S (14) 
fe=i 

where 

2711 J i-Q 

Li+L2-i l-Lmj^ 

and the integral is along the boundary, consisting of the arcs Li, L2, and Lm^- The 
parameter Afc in eq. (I14p is an integer determined by the following conditions: 



^ Re i!i£M + Afc = 0, if T Re i^l^M is an integer, (16) 
2TTi Ztti 

-1 < =F Re ^'^^^'^^'^ + Afe < 1, other cases, (17) 
2m 

where the upper sign "— " has to be taken for the starting point of a certain arc, Lj, the lower 
"+" for the end point. The condition (|17p becomes 

_l<^Rei^i^ + A,<0 (18) 



when the fundamental solution is chosen to satisfy eq. (jlip . 

By eq. ()13p . we can see that the jump G {Q) of the fundamental solution ip [Q at the two 
sides of the boundary is equal to the jump of ^' (C)j a known function defined by eqs. ([7]) and 
(l8|), along the boundary. This allows us to determine the jump G {Q by ^ {Q and then to 
determine the fundamental solution 93 {Q) by G {Q). 



2.3 Zeros 

The fugacity z is one of the zeros of the function ^ [Q) on the real axis. In principle, after 
obtaining ^ {Q), except the zero z, if we know all other — 1 zeros of ^ (C), we can write 
down the explicit expression for z as in eq. (|10p . However, the difficulty of finding the zeros 
is often the same as the difficulty of finding the zero z. That is to say, it is actually impossible 
to solve z by first solving the n^; — 1 zeros Ui. 

Alternatively, for finding the zero z, we note that eq. ([9]) is essentially an equation of z, Wj, 
and a parameter rj. Based on eq. Q, we have two possible ways to construct a set of equations 
for z, Wj, and rj: (1) Different values of C give different equations of z, Wj, and t/, and n^i + I 
different C's give a set of n^^ + 1 equations. Then, solving such a set of equations gives the 
fugacity z. (2) Deriving both sides of eq. ([9]) v times, 

, z^ = I,-- - ,n^, (19) 

C=Co 

gives equations of z, uji, and rj; for a given C, = Co) where Co is an arbitrary analytic point, 
together with eq. ([9]), we also obtain a set of + 1 equations for z, uji, and r/, whose solution 



(F^ (C) 



C=Co 



(c-^) n (c-^o^(c) 



i=l 
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gives the fugacity z. In our case, the derivations of ^ (C) and $ (C) at C = are relatively easy 
to calculate, so we will construct the equations for zeros by the second approach. 

For finding the zeros, we need to construct a set of + 1 equations. To achieve this, we 
need first to know the value of n^j, the number of zeros of ^ (C)- In our case, as the function 
^ (C) has no singularities besides the singularities on the boundary, the number of the zeros 

can be determined with the help of the argument principle: Along a contour surrounding 
the complex plane except the boundary, the change of the argument of ^ (C) is proportional 
to the number of zeros. 



3 The explicit expression for the fugacity: isochoric processes 

In this section, we will solve the explicit expression for the fugacity for hard-sphere Bose and 
Fermi gases from their equations of state directly by the method described in the above section. 
In the case of Bose gases, we will not consider the problem of phase transition. 



3.1 Hard-sphere Bose gases 

The equation of state for three-dimensional hard-sphere Bose gases, up to first order of a/ A, is 
[2] 



1 



1 



a 



^ In H (T, y ) = (2j + 1) ^ ^55/2 (z) - 2 ( j + 1) -gii^ (z) 

1 ^ 1 r a 

n = -z—lnE{T,z,V) = {2j + 1) ^ [53/2 {z) - 4 (j + 1) ^gy^ {z) 53/2 {z) 



A3 



A^ 



(20) 
(21) 



where o is the scattering length, A = h/\/2'KmkT the mean thermal wavelength, j the spin of 
the particle, and 

Qa {z) = r-— / -dt 

T{a) Jo z le* - 1 

the Bose-Einstein integral. 

In an isochoric process, the fugacity z = z {n,T) can be solved from eq. (j2ip . That is to 
say, as discussed in the above section, the fugacity z is a zero of the real function ^'i (z) = 
(l/nV) z (d/dz) InH — 1. Analytically continuing ^'i (z) to the entire complex plane gives 



(C) = (2j + 1) 



n 



A3 L' 



53/2 (C) - 4 (j + 1) -51/2 (C)53/2 (C) 



1, 



(22) 



where (C) is an analytic continuation of the Bose-Einstein integral g^ (z). Generally speaking, 
^'i (C) has more than one zero in the complex (^-plane, and the fugacity is one of these zeros 
on the real axis. 

By expressing ^1 (C) as 



n^ — l 



*i(C) = ^yi(C-^) n (C-^^)^'i(C), 



(23) 



i=l 



where z is the zero corresponding to the fugacity, (i = 1, • • • ,nuj — l) are other zeros of ^'i {() 
besides z, and r/i is a constant, we define a function <I>i ((") with no zeros. To determine $1 (("), 
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we need to analyze the behavior of the singularity of ^'i {() in the (^-plane. From eq. ()22p . we can 
see that the singularity of ^'i (C) is determined by the singularity of the analytically continued 
Bose-Einstein integrals gi/2 (C) and ^3/2 (C)- The analytically continued Bose-Einstein integral 
(z) is just the polylogarithm function, or the Jonquiere function, {(), a special case of 
the Lerch function, which is analytic in the region with the boundary along the positive real 
axis from 1 to 00 [6]. In other words, the analytically continued Bose-Einstein integral (z) 
has no isolated singularities, and all the singularities lie on the line from 1 to 00 on the real 
axis, i.e., the boundary is [1,cxd). Consequently, according to eqs. ([22]) and (f23]) . ^'i {() and 
<I>i ((") also have no isolated singularities, and all their singularities lie on the line L from 1 
to 00 on the real axis, i.e., in this case, the boundary of the analytic region consists of only 
one line. The function <I>i (C) is analytic in the region with the boundary L and everywhere 
different from zero. 

Calculating $1 {Q needs the solution of the homogeneous Riemann-Hilbert problem. Eq. 
()12p gives the relation between $1 and the fundamental solution of the homogeneous 
Riemann-Hilbert problem ipi (C). As we have chosen the fundamental solution satisfying eq. 
(jlip . the parameter /J^ in eq. (|12p can be determined by the demand that the choice of Pk 
must ensure that the behaviors of both sides of eq. (fT2l) are the same at each end. In our 
problem, the boundary of the analytic region of ^'i (C) has only one end C = 1 besides the end 
at infinity. Near the end C = 1) 

1 

i.e., ^1 (C) is divergent near C = 1 with a degree less than 1. This means /3 = 0. The 
singularities of $1 (C) and ^'i {() are the same, so <I>i (C) itself is just the fundamental solution, 
i.e., $1 (0 = 9^1 (0- 

For solving the fundamental solution ipi ((), we first need to calculate the jump, 

Gi [x) = _ = X £ [1, 00) , (25) 

(x) (x) 



*i(C)lc^i~^^, (24) 



on the boundary, according to eq. ()T3l) . 

The jump of ^'i (C) on the boundary, according to eq. ()22|) . is determined by the jump 
of the analytically continued Bose-Einstein integral gcr (C), the polylogarithm function Li^ (C)- 
The imaginary part of the polylogarithm function Li„ has a discontinuity on the boundary 

[1,00) [7J: 

ImLv {x + i6) = {InxY'-^ , 

^7 xG[l,oo), (26) 

lmLia{x - i6) = -—— {In x)" , 
T{a) 

where (5 is a small positive quantity. Therefore, the values of the analytically continued Bose- 
Einstein integral at two sides of the boundary are 

gt (x) = Lit (x) = 9a (x) ± (In x)"-^ , x E [1, 00) , (27) 

where 
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denotes the Cauchy principal value of the analytically continued Bose-Einstein integral ga (C) 
at the point x on the boundary Note that the Bose-Einstein integral go (x), i.e., the case 
of cr = 0, is an exception. gQ {() = (/ {1 — () has only one singularity ^ = 1, but has no 
singularities on the region (1,cxd) on the real axis. Accordingly, by eq. (j22p . the values of 
^'i (C) on both sides of the boundary [1, cxd) are 



(x) = (2j + 1) ^ {03/2 (x) - 4 (i + 1) ^ [gi/2 (x) 53/2 (x) - 27r] } - 1 
± i (2j + 1) -L2^Vhr^ 1 1 - 2 (j + 1) ^ ^^-03/2 (x) + 2si/2 (x) I 



(29) 



Clearly, 'iff (x) and (x) are complex conjugate to each other, i.e., ^'^ (x) = ['f^ (x)]*, and 
the jump on the boundary [l,oo) is then 

Gi (x) = exp [i2 arg (x)] , 

where the argument of "^f (x) is 

^ , , , Re (x) , , 

arg (x) = arccot h-r^- (30) 

^ ^ ^ ^ Im^+(x) ^ ^ 

Now, the fundamental solution ipi (C) can be solved by the use of eq. (jl4p directly. In our 
problem, the boundary has only one end different from infinity, C = 1) which means 

(^i(C) = e^iK)(C-l)"S (31) 

where 

,,,c, = ir..?2itM (32) 

vrA x-C 

and Ai is determined by eq. ([18]). Choosing arg^I'li'" {00) = 0, we have arg'^f (1) = —3n/2. As 
the fundamental solution satisfies eq. (jlip . the relation eq. (jl8p becomes 

- 1< -iarg^+(l) +Ai < 0, (33) 
vr 

and gives Ai = —2. Therefore, 

e7i(C) 

$1 (0 = ^1(0 = ^^^. (34) 

By eq. ()23p . we can construct a set of equations for zeros. For this purpose, we need 
to determine the number of the zeros of ^'i (C). The fact that the function ^'i ((") has no 
singularities besides the boundary [1, 00) allows us to use the argument principle to determine 
the number of its zeros in the complex (^-plane directly. Applying the argument principle along 
the contour surrounding the complex plane except the boundary [l,oo) shows that ^'i (C) has 
two zeros, denoted as z (the fugacity) and lj. Substituting eq. into eq. ([23]) with n^^ = 2 
gives 

*,(C)=,,K^ii(^e-(« (35) 
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From eq. (135]) . we can construct a set of equations for z, uj, and rji. Since in this case, the 
values of the functions in eq. (f35]) are relatively easy to be carried out at C = 0, we adopt 
the second method introduced in Sec. 12.31 Derive both sides of eq. (j35p to construct various 
equations for zeros. In this case, we need three equations for determining z, lo, and tji. When 
C = 0, we have 



571(0) 

=71 (0) 



r]izuje'"-"' = (0) , 



rjiZLue 
rjizuje'^^^^^ < 



7[{0)---- + 2 =^[(0), 

Z LO 

1 1 

- + - 

Z UJ 



j[ (Of - 2 [7; (0) + 2] 



where 



7i (0) = — 

vr ./I 



+H (o)+7;'(o) + — + 6| = ^;'(o), 

ZUJ J 

arg {x 



„n+l 



-dx. 



(36) 
(37) 



Consequently, the fugacity z can be obtained by solving eq. (|36]) : 
z = 2 { (2i + 1) + 7^0) + 2 + j {2j + if — ^ 

- (2j + 1 



{nX^y 
a' 



^ [27I (0) + 4 - ^ + 16 U + 1) ^1 + 2j'l (0) - 7; (0)' - 47; (0) ' 



-1 



(38) 



3.2 Hard-sphere Fermi gases 



The equation of state for three-dimensional hard-sphere Fermi gases, up to first order of a/X, 

ism 



p 

kf 

n 



(2j + 1) 



A3 
1 

A3 



(2j + 1)t3 h/2{z)-2j-fl^{z) 



h/2 {z) - 4i-/i/2 {z) h/2 {z) 



(39) 
(40) 



where the Fermi-Dirac integral 



fAz) 



Via) 



4.<T-l 



z-^e* + 1 



-dt. 



In an isochoric process, the fugacity z = z {n,T) can be solved from eq. (j40p . The fugacity 
z is a zero of the complex function 



^2 (C) = (2j + 1) 



n 



A3 



/3/2(C)-4j-/i/2 (0/3/2(0 



1 



(41) 



on the real axis, where (0 is the analytically continued Fermi-Dirac integral. The singularity 
of ^'2 (C) is determined by the singularity of the analytically continued Fermi-Dirac integral. 
The analytically continued Fermi-Dirac integral here is the polylogarithm function — Lig- (— C); 
which is analytic in the region with the boundary [—1, —00). 
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By the argument principle, along the contour surrounding the complex plane except the 
boundary [— 1,— oo), we can determine that ^2 (C) has only one zero on the real axis, the 
fugacity z. Accordingly, ^'2 (C) can be expressed as 

*2 (C) = r?2 (C - ^) $2 (C) , (42) 

where the function ^*2 (C) has the same singularities as those of ^'2(C)) lying on the line 
[— 1, —00), and has no zeros. 

As stated above, in our case, the fundamental solution is chosen to satisfy eq. (Ilip . Near 
the end different from infinity of the boundary [—1, —00), C, = —1, we have 



'I'2 (C)lc^-l ~ ^ 



1 



+ 1 



(43) 



The degree of divergence of ^'2 (C) a-t the end C, = —\ being less than 1 implies that the 
function $2 (C) itself is a fundamental solution of the homogeneous Riemann-Hilbert problem, 
i.e., $2 (C) = ^2 (C)) where ^f>2 (C) denotes the fundamental solution. 

To solve the fundamental solution ip^ (C)) we need to know the jump on the boundary. From 
eq. ([I3D, we have 

- Hi -IM' 

and from eq. (j26p . we have 

ft (x) = -Lit i-^) = U (x) T ijr^ [In i-x)r' , (45) 

where 

U ix) = :p^P / _i , ^ dt, a + (46) 

is the Cauchy principal value of the analytically continued Fermi-Dirac integral. Notice that 
X < 0. Then, the value of ^'2 (C) at the two sides of the boundary can be obtained directly: 

^± (x) = (2j + 1) {f3/2 ix) - 4j^ [fi/2 (x) f3/2 (x) - 27r] } - 1 



(47) 



(x) and ^'2^ (x) are complex conjugate to each other, so the jump of ^'2 (x) on the boundary 
[—1, —00) is 

G2 {x) = exp [i2 arg (x)] , 

where the argument arg (x) = arccot [Re (x) / Im ^'^ (x)] . Noting that in this case, 
there is only one end different from infinity of the boundary, C, = —1, we can write down the 
fundamental solution from eq. (jl4p : 



<^2(C) = e^^(^)(C + l)"% (48) 

where 

T^J-l X-C 
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and A2 will be determined by the condition ()18p . Choosing arg (—00) = gives arg ^'^ (—1) = 
— 7r/2. As the fundamental solution satisfies eq. ([TT]) . the condition (fTSll becomes 

-1< --arg^'+(-l) + A2 <0, (50) 

TT 

and gives A2 = — 1- Then, 

e72(C) 

'52(C)=9'2(C) = ^. (51) 

Substituting eq. ([5T]) into eq. (|12]) . we have 

*2(C) = %^e^^(^). (52) 

For determining the zero z and the constant 772 , we need two equations. By setting C = 0> sq. 
(j52|) and the derivative of eq. ([52|) give these two equations. Consequently, we have 

(53) 

4 The explicit expression for the fugacity: isobaric processes 

In the above section, we calculate the explicit expression for the fugacity in isochoric processes. 
In this section, we will calculate the explicit expression for the fugacity for hard-sphere quantum 
gases in isobaric processes. Isobaric processes are of great significance to the problem of phase 
transition. In an isobaric process, the fugacity is determined by eq. ([T]). 



(2j + 1) 



n 



A3 



72(0)-l 



4.1 Hard-sphere Bose gases 



The equation of state for hard-sphere Bose gases is given by eqs. ()20p and (|2ip. In isobaric 
processes, from eq. ([20]) . the fugacity z is a zero of the complex function 



(C) = (2i + 1) 



kT 

PA3 L- 



55/2 (C)-2(i + l)-<73'/2 (0 



1 



(54) 



on the real axis. The boundary of the analytic region of ^'5 (C) is the line [1,cxd), and ^'5 (C) 
can be written as 

^5(C)=%(C-^)$5(C). (55) 

$5 ((") has no zeros and has the same singularities as ^5 ((")• Since ^'5 ((^) is convergent at the 
unique end different from infinity, (" = 1, of the boundary, $5 (C) = '^^ (C) is just a fundamental 
solution. 

At the two sides of the boundary [l,oo), by eq. (p7|) . we have 



kT 



(x) = (2j + {g5/2 (x) -2{j + l)j 



03/2 (x) 



47r Inx 



1 



± i {2j + 1) -—A^Vln 



PA3 



^ Inx - 2 (j + 1)^03/2 (x) 



(56) 
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The jump is then 



G5 (x) = = = exp [i2 arg ^+ (x)] , x G [1, 00) 



(57) 



(x) (x) 

where arg (x) = arccot [Re (x) / Im (x)] . The fundamental solution is consequently 

c^5(C) = e^^K)(C-l)'% (58) 

where 



x-C 

The fundamental solution is chosen to satisfy eq. ()lip . We choose arg (00) = 0, then 
arg^f^ (1) = —IT, thus the condition 



gives A5 = —1, and so 



--Rearg^'+(1) + A5 = 

TT 



^5 (C) = %^e^^(^) 



(60) 



(61) 



C-1 

By setting C = 0, eq. (f6T|) and its derivative give two equations to determine the zero z 
and the constant r/5, and the solution is 

1 -1 



(2j + l);^+75'(0) + l 



(62) 



4.2 Hard-sphere Fermi gases 

Similarly, for the case of hard-sphere Fermi gases, according to eq. ([39|) . we define 



*6 (C) = (2j + 1) 



kT 
PA3 



h/2 (C)-2j^/3% (C) 



1, 



(63) 



which has a unique zero on the real axis, the fugacity z, and the boundary of its analytic region 
is [—1, —00). ^6 (C) can be expressed as 



*6 (C) = % (C - Z) 1>6 (C) • 

As (C) converges at the unique end different from infinity of the boundary, C = 
(pQ (C) is a fundamental solution. 

At the two sides of the boundary [—1, —00), 

kT ( a r i ~i 

^± (x) = (2j + 1) {f5/2 {x) - 2j- [f2/2 (x) - 47r In (-x)J | - 1 



(64) 

-1,$6 (0 = 



T iiV^ (2j + 1) Vln(-x) 



^ln(-x) - 2j^f3/2 (x) 



and the jump is 



Ge (x) = 4^ = ^^[4 = [^2arg (x)] , x G [-1, -00) , 

IPQ (x) (x) 



(65) 



(66) 
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where arg^'^j" (x) = arccot [Re^^ (x) /Im\I'^ (x)] . Thus, the fundamental solution is 

(^6(C) = e^«(^)(C + l)'% (67) 

where 

The fundamental solution is chosen to satisfy eq. Choosing arg ^'g' (— cxd) = gives 

arg^g" (—1) = —IT. Then, the condition (fTUj) 



-Rearg^+(-l) + A6 = 



TT 



gives Ae = —1- Consequently, 



^ ~ ^ .76(0 



*6 iO = V6-^e 



(69) 



(70) 



By setting = 0, eq. ([70]) and its derivative give two equations for z and the constant 
and the solution of the fugacity is 

-1 



kT 

(2j + l)pi^+76(0)-l 



(71) 



5 Conclusions 



To sum up, in this paper, we present the explicit expression of the fugacity for three-dimensional 
hard-sphere Bose and Fermi gases in isochoric and isobaric processes. The method is to convert 
the problem of solving the fugacity from the equation of state into the problem of finding the 
zero of a complex function based on the homogeneous Riemann-Hilbert problem. This method 
is introduced by Leonard for treating ideal quantum gases f5]. Concretely, in this treatment, 
one can solve the fugacity z from the relation 

^(C)=r?(C-z) n (C-^0 „. — , 

fc=i 

where z and uji, i = 1, • • • , tIi^ — 1, are the zeros of the complex function ^' (C) which is 
constructed from the equation of state of quantum gases, is the fc-th end that is different 
from infinity of the boundary of the analytic region of ^ (C), n^^ is the number of the zeros, 
n is the number of the ends different from infinity of the boundary, ip is a fundamental 
solution of the homogeneous Riemann-Hilbert problem, and fi^ is a constant chosen to ensure 
that if (C) is a fundamental solution and is determined by the behavior of ^ ((^) near the end 
Cfc. By deriving both sides of this equation at a given point or setting various values of C, one 
can obtain a set of equations of z, Wj, and the constant r]. Solving such a set of equations gives 
the fugacity z. 

The key steps in this treatment are (1) analytically continuing the real function ^ (z) whose 
zero corresponding to the fugacity z to the whole complex plane, which gives the complex func- 
tion ^ ((^), (2) finding the fundamental solution, and (3) determining the number of zeros of 
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^'(C)- At the first step, in the present case, the analytic continuation of ^ (z) rehes on the 
analytic continuation of Bose-Einstein and Fermi-Dirac integrals; the analytically continued 
Bose-Einstein and Fermi-Dirac integrals are the polylogarithm (Jonquiere) functions, Liu (C) 
and —Liu ("Oj respectively. At the second step, based on the result of the analytic continua- 
tion, we can analyze the singularity structure of ^ determine the boundary of the analytic 
region of ^ (C), calculate the jump of ^ at the boundary, and then solve the fundamental 
solution <^{C)- At the third step, though in general by the argument principle, one can only 
determine the difference between the number of zeros and the number of isolated singularities, 
one can determine the number of the zeros by use of the argument principle due to the fact 
that ^ {() has no isolated singularities in the present case. 

It is worthy to point out that the explicit expression of the fugacity for the Bose case given 
in the present paper, especially the result for isobaric processes, can be directly applied to 
analyze the phase transition of the hard-sphere Bose gas system, which is an important problem 
in recent times [9]. The reason is that the fugacity, and therefore the chemical potential, 
plays a central role in the theory of phase transitions, according to Erenfest's theory of phase 
transitions. In the Erenfest classification, the order of a phase transition is defined by the order 
of the discontinuities in the derivatives of the Gibbs free energy, G = Nfi. A detailed discussion 
on this problem will be given elsewhere |10] . 

We are very indebted to Dr. G. Zeitrauman for his encouragement. This work is supported 
in part by NSF of China under Grant No. 10605013 and the Hi- Tech Research and Development 
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